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' Abstract. We discuss some basic issues that arise when one attempts to model 

quantum mechanical systems on a computer, and we describe the mathematical 



structure of the resulting discretized cannonical commutation relations. The C*- 
algebras associated with the discretized CCRs are the non-commutative spheres of 
Bratteli, Elliott, Evans and Kishimoto. 
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1. Introduction. 

We discuss some basic issues associated with the problem of modeUing quantum 
mechanical phenomena on a computer. Before setting out to write code for any 
program of this kind, one first has to recast the quantum mechanical Hamiltonian 
into a form appropriate for doing numerical calculations; at the same time, one 
has to retain the essential features of quantum mechanics (i.e., the uncertainty 
principle). Reflection on this problem shows that one must confront a more basic 
issue, namely that of "discretizing" the canonical operators P, Q which satisfy the 
commutation relations 

PQ-QP=\l. 

% 

For simplicity, we consider one-dimensional quantum systems. We will show 
that there is a natural (perhaps one could say an inevitable) way to discretize the 
canonical operators (P, Q) so that (a) the discretized pair conforms to the basic 
principles of numerical analysis, and (b) the uncertainty principle is preserved. 
We will see that there is a one-parameter family {Pa,Qa)i a > 0, of discretized 
canonical operators. The parameter a is the square of the numerical step size. For 
each a, Pq, and Qa are bounded self-adjoint operators which obey an altered form 
of the canonical commutation relations. 

If a is not a rational multiple of 27r then the unital C*-algebra generated by 
{Pai Qa} turns out to be one of the non-commutative spheres of Bratteli et al [1] , [2] , 
and hence is a simple C*-algebra which has a unique tracial state. This C*-algebra 
contains the discretized Hamiltonian, and that opens the possibility of computing 
its spectrum. While much remains to be learned, there are examples in which the 
spectrum of the discretized Hamiltonian is totally disconnected. In particular, these 
C*-algebras contain many projections. We also show how each of these C*-algebras 
arises as the enveloping C*-algebra associated with an appropriately discretized 
version of the Weyl canonical commutation relations. Thus, it seems appropriate to 
conclude that non- commutative spheres will arise naturally in any serious attempt 
to model quantum phenomena via computer. 

It is unlikely that there is a direct connection between these non-commutative 
spheres and the quantum spheres of Woronowicz and Podlcs [6]. I should also 
point out that this work is in progress, and that full details will appear elsewhere. 
I would like to thank Paul Chernoff and Alan Weinstein for objecting strenuously 
to an earlier version of this paper; their constructive criticism has caused me to 
think more carefully about the exposition, with the result that the paper is now 
better than it was. Whatever deficiencies that remain are the sole responsibility of 
the author. 

2. Background. 

In order to tie the discussion to a specific issue, we digress briefly in order to 
describe the kind of computing problem which has motivated this work. 

Consider the one-dimensional anharmonic oscillator. Apart from a choice of 
units, the relevant equation of motion is 

2.1 x + x + gx-^^O. 

Here, is a nonnegative constant and x = x{t) is a real- valued function of time t. 

In physical terms, one has a unit mass attached to one end of a nonlinear spring. 
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length L, the restoring force is L + gL^. The spring is hnear (i.e., it obeys Hooke's 
law) when ^ = 0. 

In order to quantize this system, one passes to the Hilbert space L^(R) in the 
Schrodinger picture. Given an initial wave function / G L^(M) at time zero, basic 
principles dictate that / evolves in time according to 

f{t) = e^*^/, t e R, 

where H is the quantum mechanical Hamiltonian. In the case of the anharmonic 
oscillator, H is the operator 

where P = —id/dx^ Q is multiplication by x, and v{x) = x'^/2 + gx'^/A is the 
potential function associated with the differential equation (2.1). 

It is not hard to "draw" solutions of the classical equation (2.1) using a modern 
desktop computer. I have done this with good results for systems of up to 50 inter- 
acting anharmonic oscillators. For example, even with small systems (four masses, 
say) one can observe interesting chaotic behavior resulting from the nonlinearity 
of the system for > 0. On the other hand, it is less clear how one might use a 
computer to create an interesting graphical representation of the quantized system. 
What I have in mind is to draw pictures of a quantized system of interacting an- 
harmonic oscillators which resemble the pictures described above, at least insofar 
as that is possible. 

With every quantum system there is an associated stationary stochastic pro- 
cess whose state space is the configuration space of the classical system ([4], [5]). 
Without going into the details of how this process is constructed (see section 5), 
let me just say that my intention is to construct "typical" sample paths for this 
quantum process and draw them on a computer screen. While I am not yet finished 
writing the code for this program, I am nevertheless convinced by the results below 
that it is possible to do it. Practical obstructions, relating to the runtime for such 
programs, remain to be seen. 

3. Discretizing P and Q. 

We turn now to the specific issues that arise when one seeks to build numerical 
approximations to a quantum system. One must first replace the Hamiltonian 

H^^P' + viQ) 

with an approximation to it which is appropriate for doing numerical studies. Let 
us first consider the momentum operator P = —id/dx. Letting V be the one- 
parameter unitary group of translations 

Vtf{x) = f{x-t), 

then of course P is the infinitesimal generator of V 



Equivalently, we have 



P = liTTI -(V, - 1) 
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in an appropriate sense on the domain of P. Conventional wisdom indicates that 
P should be replaced with a difference operator, the simplest one being 

Pr = -{Vr-l). 
IT 

Here, r is the numerical step size. While this operator is certainly bounded, it is the 
wrong choice for two reasons. It is not self adjoint, and for small r it approximates 
P only to terms of order 0(t), which is not very good. Thus it is better to use the 
real part of the above operator, namely 

This is a bounded self adjoint operator which approximates P (on smooth test 
functions) to terms of order O(t^). Using the fact that Vr = e^'^^ we have the 
equivalent formula 

3.1 Pr = -sin(TP). 

T 

Considering that the function 

sin(ra;) 



/(^) = 

•/ 

approximates the function g[x) — x very closely in a fixed interval —M < x < +M 
when T is small, one can expect good numerical results from P^- for small values of 

T. 

We must now approximate Q in a similar fashion, but in doing so we must be 
careful to preserve the essential features of quantum mechanics. That means the 
uncertainty principle. In fact, the uncertainty principle for the approximating pair 
{Pr, Qt) is achieved by choosing 

Qr = F-^PrF, 

where F is the Fourier transform operator on (R) . In order to discuss our reasons 
for making this choice, we digress momentarily to a more coordinate-free setting. 

Suppose that one has a pair of one-parameter unitary groups U, V which act on 
the same Hilbert space H in such a way that VsUt — e^^^UtVs, and such that the 
set of operators {Ut} U {Vs} is irreducible. Then the self-adjoint generators P, Q 
of these two groups 

Vt = e**^, Ut = e^*^ 

obey the commutation relation cited in the introduction. Choose a unit vector / 
in the underlying Hilbert space H. Then the probability measure associated with 
the observable Q and the state / is the probability measure defined (on the Borel 
subsets of M) by 

q{S)^<EQ{S)fJ>, SCR 

where Eq is the spectral measure of Q. Similarly, the probability distribution 
associated with P and / is 



NON-COMMUTATIVE SPHERES AND NUMERICAL QUANTUM MECHANICS 5 



The uncertainty principle for the pair of observables {P, Q) makes the following 
qualitative assertion: if a wave function / represents a physical state such that the 
measure q appears highly concentrated near a point, then the measure p must ap- 
pear "spread out" . A similar assertion applies, of course, with p and q interchanged. 

It is not hard to reformulate this somewhat vague assertion entirely in terms 
of the Fourier transform operator. To see this, notice first that there is a unitary 
operator F e B{H) such that 

3.2 P = P-^QF, -Q = P-^PP, 

Indeed, since the pair (P', Q') defined by P' = — Q, Q' = P is irreducible and 
satisfies the same commutation relation ( equivalent ly, U-gVt = e^^^VfU-s, s,t G M), 
the Stone- von Neumann theorem implies that the two pairs (P, Q) and (P', Q') are 
unitarily equivalent, and P may be chosen as a unitary operator implementing 
this equivalence. Moreover, because of the irreducibility hypothesis P is uniquely 
defined by (3.2) up to a scalar multiple of absolute value 1. It follows from the 
relation P = P~^QP and elementary functional calculus that 

Ep{S) = P-^Eq{S)P, SCR, 

and hence the formula for the measures p and q can be rewritten in terms of P as 
follows 

qiS) = <EQiS)fJ> 

p{S) = <EQ{S)Pf,Pf>, SCR. 

Applying the Stone-von Neumann theorem once again, we may assume (after 
replacing {P,Q) with a unitarily equivalent pair) that H — L^(M), P — —id/dx, 
Q = multiplication by x, and that / is a normalized function in L^(R). Notice too 
that, after this change of coordinates, P must be a scalar times the usual Fourier 
transform operator. Indeed, for this choice of P and Q acting on L^(]R), the Fourier 
transform operator Pq is well-known to satisfy (3.2), i.e., 

PoP = QPo, PoQ = -PFo. 

Hence by irreducibility we must have P = XPq where A is a complex number of 
absolute value 1. Thus in these coordinates the preceding formulas for p and q 
become 

q{S)^ I \f{x)\''dx 
Js 

p{S) = [ \f{x)\^dx, SCR 
Js 

f denoting the Fourier transform of the function /. 

At this point, one can recognize the uncertainty principle as reflecting a familiar 
property of the Fourier transform. Indeed, if the measure q is highly concentrated 
near a point xq then so is the wave function /; hence its Fourier transform / 
is highly spread out (approximating a constant multiple of the pure exponential 
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More generally, if one starts with any bounded self-adjoint operator Pi in the von 
Neumann algebra generated by all translation operators on (R) , then one can ex- 
pect some form of the uncertainty principle for the pair of operators {Pi,FPiF-^), 
at least for bounded operators Pi which are reasonable approximations to P. In 
more detail, we first find a bounded real valued Borel measurable function (/) such 
that Pi = (j){P). Standard functional calculus implies that Qi = FPiF~^ = 
F(f){P)F~^ = (j){Q). Thus if / is any normalized function in L^(R) then we have 



for any Borel set in M, (j)~^S denoting the inverse image of the set S under 0. It 
follows that the probability measures associated with the observables Qi, Pi and 
the state defined by / are given by 



In the most degenerate cases where (p is nearly constant, these formulas have 
little content. However, if — x\ is small for all x in a very large interval, 

then these two measures closely approximate the measures p, q of the preceding 
paragraphs, and hence they can be said to obey the uncertainty principle. 

Returning now to the problem of approximating Q, wc conclude from these ob- 
servations that in order to preserve the uncertainty principle for the approximating 
operators (Pt-,<5t), one should set 



where F is the usual Fourier transform. Notice that there is no arbitrariness in 
the choice of Qr', once we have settled on a discretized approximant Pr to P, Qr 
is forced upon us by requiring that the uncertainty principle should be preserved. 
Equivalently, the definition of Qr is 



< EQ^{S)fJ) = < EQ{cp-'S)fJ>,and 

< Ep,{S)fJ) = < Ep{4>-'S)f,f>=< EQ{r'S)Ff,Ff > 




Qr = FPrF 



--1 



3.3 



Qr = 



sin(r(5) 



T 



In terms of the one-parameter group of multiplication operators 



UJ{x) = e^'-fix), 



this definition can also be written 



Qr^^{Ur-U.r), 



T being the same step size used for Pr- 

Thus for every r > 0, we have a pair of bounded self adjoint operators {Pr, Qr)- 

\-HT„ ;ii r^* I Ti r\ \ c — r^* „ j- — ] u,, o i r\ 



NON-COMMUTATIVE SPHERES AND NUMERICAL QUANTUM MECHANICS 7 



Finally, our approximation to the original Hamiltonian is defined as 

We abuse notation slightly here, but since we will have no need to refer back to 
the original unbounded Hamiltonian it will cause no problem. H is a, bounded 
self adjoint operator for every continuous real-valued function v, and it obviously 
belongs to C*{Pr,QT)- 

defines a uniformly continuous one-parameter group of inner automorphisms of this 
C*-algebra. 

Let us now recall the definition of non-commutative spheres Ba [1], [2]. Let 
q; e M and let U,V he a pair of unitary operators which satisfy 

VU = e"^UV. 

The C* -algebra Aa generated by U and V is well known to be independent of the 
particular choice of U and V, is simple, and has a unique tracial state. 

Noting that the pair {U~^,V~^) satisfy the same commutation relations as 
{U,V), we can define a unique *-automorphism a of Aa by requiring 



3.4 a{U) = U-\ a{V) = V'K 

We have cr^ = id, and hence cr defines a Z2-action on Aa- When a is not a 
rational multiple of tt, Ba is defined as the fixed subalgebra 

Ba = {a e Aa ■ cr(a) = a}. 

Theorem 3.5. For every positive r such that t^/tt is irrational, C*{Pr,QT) is 
isomorphic to B^-^ . 

Remarks. In case ck/tt is irrational, it is well known that the unique trace on Aa 
gives rise to a representation of Aa which generates the hyperfinitc Hi factor R, 
and one can show that the weak closure of the d-fixed subalgebra is a subfactor 
of R having Jones index 2. Recalling that such subfactors of R are unique up to 
conjugacy (see [3], pp. 1-2 for more detail), we can recognize this as a very stable 
invariant for the embedding of C* {Pr, Qt) in at the level of von Neumann 
algebras. Thus, we may think of C*{Pt-, Qt) as being isomorphic to an "index two" 
C*-subalgebra of the irrational rotation algebra Ar^ ■ 

We also remark that since the step size r must be chosen to be a small rational 
number in any realistic computational setting (for example, r = 10~^ is typical 
for me), the C*-algebras C*(Pt-,Qt-) that actually arise will always satisfy the 
hypothesis of Theorem 3.5. 

Let tr be the tracial state of C*(Pt-, Qr)- Then for every (3 > there is a KMS 
state ojj3 associated with the group 7 and the value up is defined on C*{Pt-, Qt) 

by 

tr{e-f^^a) 
'"^^"^ = tr{e-P^) ■ 

The GNS construction applied to ujjj leads to a representation of C*{Pr, Qt) which 
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4. The discretized CCR algebra.. 

We now give an alternate description of C*{Pri Qr) as the universal C*-algebra 
associated with a discretized version of Weyl's form of the canonical commutation 
relations. Significantly, it is this realization of C*(Pr,<5r) that is most useful for 
doing numerical calculations. Fix o: > 0, not a rational multiple of n. Consider 
the discrete abelian group G = Z x Z, and let a; : G x G — > T be the bi-character 
defined by 

a;((m, n), (p, q)) = e*("f-^«)«/2. 

By a representation of the discretized canonical commutation relations we mean a 
uniformly bounded family {Dx '■ x e G} of bounded self adjoint operators on a 
Hilbert space which satisfy 

4.1 DxDy = u){x, y)Dx+y + uj{y, x)Dx-y. 

Remarks. Formula 4.1 is a generalization of the elementary trigonometric identity 

2 cos AcosB = cos{A + B) + cos{A - B), 

in which phase shifts have been added via the cocycle u. For example, the real- 
valued function Z> : G — > M defined (for a and fixed) by 

D{m, n) = 2 cos(mcr + ncf)) 

satisfies 4.1 for the trivial cocycle u = 1. It is also clear that 4.1 is closely related 
to Weyl's formulation of the canonical commutation relations. Finally, it can be 
shown that any bounded family of self-adjoint operators satisfying 4.1 must in fact 
obey \\Dx\\ < 2; thus there is a natural construction which produces a universal 
enveloping C* -algebra C*{G,a) for the relations 4-i- 

It can also be shown that the original operators Pr, Qt serve to define a unique 
*-representation tt of C*{G,t^). it is specified uniquely by requiring that 

4.2 7r(i:»(l,0)) = 2tP^, 7r(i:>(0, 1)) = 2tQ^. 

While a fairly explicit formula can be given for each of the operators n{D{m, n)) in 
terms of Pr and Qr, we shall not do so here. Now we are in position to formulate 
the second characterization of C*{Pr, Qt)- 

Theorem 4.3. The representation n defined by 4-2 gives rise to an isomorphism 
of C* -algebras 

C*(Z X Z,T^) ^ C*{Pr,Qr)- 

5. The canonical stochastic process. 

We now indicate how one can construct a stochastic process using the discretized 
operators Pr, Qr from section 3. Let v be any continuous real- valued potential 
function, and set 

H=^P^ + viQr). 
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Theorem 5.1. Let fo,fi,---,fn be nonnegative continuous functions of a real 
variable, and let ti,t2, ■ ■ ■ ,tn > 0. Let tr be the tracial state on C*{Pr, Qt)- Then 

tr{fo{Qr)e-''''fl{Qr)e-'''' ■ ■ ■ e-'-" fn{Qr)) > 0. 

Using Theorem 5.1, we can make use of apparatus developed by Klein and Lan- 
dau [4] as follows. Let a >0 and let tr be the tracial state of Theorem 5.1. 

Theorem 5.2. Fix /? > 0. There is a real-valued stochastic process {Xt : < t < 
P} which is defined by the conditions 

E{f,{Xt,)f2iXt,)...fn{XtJ) = 

tr{e-PH) 

for every < ti < t2 < ■ ■ ■ < tn = every fi,f2,---,fn e C'oW, and every 
n= 1,2,.... 

Remarks. This process {X^ : < t < P} is stochastically continuous and periodic 
with period /? in that Xi3{u)) = Xq{u!) almost surely. Moreover, if one extends it 
in the obvious way to a periodic process (with period f3) defined for all t G M, 
the resulting process is stationary and obeys a form of reflection-positivity (i.e., 
Osterwalder-Schrader positivity) which is appropriate for periodic processes. 
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